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1. Introduction 

Self- avoiding walk models are of interest as models of polymers in solution pQ as well as 
realisations of 0(n) spin models in the limit n —>■ 0[2 . 

The self-avoiding walk problem is a difficult problem numerically because of the 
intrinsic long ranged nature of the problem; self-avoiding walks are extended objects. 
The numerical methods most used to date for the study of such models are Monte Carlo 
simulations|3], series expansions 4 and transfer matrix calculations El 13 El El CHI- 
Monte Carlo simulations are hampered by problems of trapping in the low 
temperature, or equivalently large interaction, phases. These may be relieved using 
the fluctuating bond method [HI E], which also enables the investigation of the dense 
walk phases. This method has the drawback of modifying the model in a way which is 
not always convenient, particularly if frustration effects are of interest. Current Monte 
Carlo methods, such as the PERM method[T2*] and parallel tempering^], relieve these 
difficulties to some extent, but at the cost of computational complexity. 

Series expansions are limited by the lengths of walks that may be enumerated. In 
practice these are limited to a few tens of lattice steps. Very technical methods may 
be used to achieve longer series expansions of the partition function (of the order of 60 
lattice steps), and in some cases the estimate of the critical point is very impressive |14j. 
If correlation length exponents are required, then the range of walk lengths open to 
enumeration are of the order of 30 steps [T^]. This becomes increasingly limiting as the 
complexity of the model increases. 

Transfer matrices have the advantage over Monte Carlo simulations of giving 
numerically exact values of the free energy, correlation lengths, polymer densities,etc, 
for strips of infinite length, but of finite width. Finite size analysis methods then allow 
the extrapolation of the finite size estimates. The widths of lattice accessible are limited: 
for the self-avoiding walk model widths up to 12 or 14 are possible (see for example |16j), 
but this again drops as the complexity of the model increases (see for example |17j). It 
often occurs that the number of data points one may use in finite size analysis becomes 
dangerously small for finite size scaling techniques to be of use, particularly if odd/even 
parity effects are present. 

In this article we apply a recent extension of the corner transfer matrix 
renormalisation group (CTMRG) method JH] to a class of self-avoiding walks introduced 
by Blote and Nienhuis, which we shall refer to as the vertex-interacting self-avoiding 
walk model, to reflect the fact that the monomer-monomer interactions are added for 
doubly visited sites (for a full description of the model see Section l2~T| . We will show 
that the CTMRG method extends the possibilities of transfer matrices to larger system 
sizes, at the cost of introducing a phase space truncation scheme. The accuracy of the 
results is however controllable through a parameter m corresponding to the number of 
states kept from one iteration to the next. 

The other important result in this paper relates to the vertex-interacting self- 
avoiding walk model. Whilst it has the standard collapse transition as a special case 
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(where the walk is obliged to turn after each step), it is believed to have a generic collapse 
transition in a different universality class^H]- We show that this collapse transition has 
a first order character, despite the fact that there exists critical exponent v = 12/23 
which has been exactly determined by Warnaar et al [19J. 

2. The corner transfer matrix renormalisation method 

Transfer matrices, as a numerical method, have a major advantage over Monte Carlo 
methods in that they sum over all possible configurations, and so are not hampered by 
trapping effects at low temperatures (or high densities in polymer problems) which 
plague Monte Carlo simulations of models with competing interactions (frustrated 
models). 

The major inconvenience is the large amount of computer memory required to 
store the matrices, limiting the maximum lattice width in practical calculations. This 
is alleviated in part by finite-size analysis methods which work well 20J, at least for the 
simpler models, coupled with efficient extrapolation methods |2*T]. 

In two-dimensional classical spin models, or one-dimensional quantum spin models, 
there exist iterative approximation schemes which, by successive pruning of phase space, 
permit calculations for much larger lattice widths for given computer resources. These 
methods are the Density Matrix Renormalisation Group method (DMRG) [22 123 and 
the Corner Transfer Matrix Renormalisation Group method (CTMRG) |24j. These 
methods have little recourse to the renormalisation group, and are better thought of 
as iterative matrix approximation schemes. The basic idea is that the (approximate) 
transfer matrix for a lattice of size 2N is calculated from the transfer matrix for a lattice 
of size 2N — 2 by inserting extra lattice sites. A change of basis is then performed, such 
that the transfer matrix may be projected onto a smaller basis, with the smallest possible 
loss of information. 

In a recent paper we showed how the CTMRG method, first introduced for classical 
spin systems by Nishino[2"^]. may be applied to self- avoiding walk type models [TH]. This 
requires a non-trivial modification of the method since the iteration scheme requires all 
interactions to be local. The self-avoiding walk problem, by its very nature, is non-local. 
In the remainder of this section the details of the method, absent in reference j!8j. are 
provided. In particular we show how this non-locality may be overcome, and describe 
the implementation of CTMRG for this class of problem. 

2.1. The self-avoiding walk, the 0(n) model and vertex models with complex weights 

It is known that spin models with 0(n) invariant spins have high temperature expansions 
consisting of closed loops. The self-avoiding loop (SAL) model is related to the high 
temperature expansion of the 0(n = 0) model[2]. In a similar way, the high temperature 
expansion of the susceptibility leads to a model with one open path, and a gas of 
loopsj22]. The n — > limit eliminates the loops to leave a self-avoiding walk. The 
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partition function for a self-avoiding walk is thus equivalent to the susceptibility of an 
0(n — 0) spin model. 

The connection between the 0(n) class of models and better known models is as 
follows: n = 1 corresponds to the Ising model, n = 2 to the XY model and n = 3 to the 
classical Heisenberg model. The partition functions of these models are given by: 

Z = E ex P ( \fi J E ■ ? i J > W 

{*} V <M> / 

where refers to a sum over nearest neighbour spins. Whilst the high-temperature 
expansion of the Ising model is relatively easy, the expansion for other values of n is a 
little cumbersome. In what follows we shall consider a slightly different 0(n) invariant 
model, due to Nienhuis ffiH l2"T] . with the partition function given by 

zo {n) = e n f 1 + k e W) > ( 2 ) 

K } V a=l J 

where the spins, sf = ±1, are placed on the lattice bonds. A diagrammatic expansion of 
Equation |21 follows if we identify the 1 as the weight of an empty bond between the sites 
% and j and the K as the weight of an occupied bond. Since (sf) = and {sf) 2 = 1, 
the only terms in the expansion of Equation |2] which survive the sum over {sf} are 
those terms consisting of closed loops where the 'colour' a is conserved around each 
loop. Different loops may have different 'colours'. The sum over a gives a factor n for 
each loop. An equivalent expression for Equation El may now be given in terms of loop 
graphs Q: 

Zo {n) = Y.^ m K^\ (3) 
Q 

where I is the number of loops and h is the number of occupied bonds. The partition 
for self- avoiding loops is then just 

Z S al = lim -Zou)- (4) 

The parameter n is now a fugacity controlling the number of loops, and in this 
representation need no longer be taken as an integer. This fugacity corresponds to a long- 
range interaction, since the loops may be of any size. This non-locality is undesirable 
for our purposes. We would like to express n as a product over local weights despite the 
variable size of the loops. The only real numbers with this property are and 1. This 
leads us to introduce complex local weights, Wi. If \wi\ = 1, the product of the Wi around 
a loop will also be of modulus 1. The final weight must, however, be real. This may 
be achieved if for every complex weight, we also have its complex conjugate. Each loop 
may be followed clockwise or anticlockwise. The loops in Equation El are not oriented, 
but may be oriented by associating 2 l ^> oriented graphs with each non-oriented graph. 
The loop fugacity is not required to be the same for the two orientations of the loop, and 
so a fugacity n + is associated with (anti) clockwise oriented loops. The partition 
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function may then be written [251 12*7] 

Zo[ n) = EK + n -W ( 5 ) 

Q 

= J2^n l ZK b (6) 
Q> 

where Q' is the set of oriented loop graphs and l + (/_) is the number of (anti) clockwise 
oriented loops. On the square lattice there must be four more corners with one 
orientation, compared to the other orientation, in order to close a loop. Using this fact 
we may set Wi = exp(i6/4) for a clockwise corner, W{ = exp(— iO/4) for an anticlockwise 
corner and Wi = 1 otherwise. This leads to n + = exp(i8) and n_ = exp(—i8), and hence 
n = 2cos(#). The oriented loop factor has now been broken up into local weights and 
the partition function may be rewritten in terms of a vertex model[2<""| I29j: 

En** w 

g> i 

where Vi is the weight of the vertex found at site i. The possible vertices are shown 
in Figure The derivation given here is only for the simplest case, but we may freely 
change the weights of the vertex configurations in order to generate different interactions 
in the original model. Notably, to have exactly the standard self-avoiding walk model it 
is necessary to set the weights of vertices 14 to 19 to zero, and set n = 0, or equivalently 
set 9 = 7r/2. 

A vertex model may be thought of as a spin model where the spins live on the bonds 
of the lattice. Each spin has three states corresponding to the empty lattice bond, a 
bond oriented in the positive x (y) direction and a bond oriented in the negative x 
(y) direction. The vertex weights may be written as a vertex function Wfa, o~j, 07., 07), 
where the tensor elements correspond to the vertex weights of Figure H The vertex 
function was defined using the following convention: 

1 if the bond has an arrow to the right or upwards, 

(j=^0 if the bond has no arrow, (8) 

-1 if the bond has an arrow to the left or down. 



2.2. The vertex-interacting self-avoiding walk model 

In this section we apply the CTMRG method explicitly to the vertex-interacting self- 
avoiding walk model introduced by Blote and Nienhuis|2""""t !3Uj . This model consists of a 
random walk which is allowed to collide at a site, but is not allowed to cross or occupy 
a bond more than once. An interaction energy e < is introduced for each collision, 
which replaces the attractive nearest-neighbour interaction present in the standard 0- 
point model. A step fugacity, K, is introduced in order to control the average length of 
the walk. A stiffness is also introduced by weighting sites sitting on a bend in the walk 
differently from sites sitting on straight portions of the walk. Following the convention 
of Blote and Nienhuis|3t)J, we choose to add a weight p to the straight segments of the 
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Figure 1. The 19 allowed vertices in the Nienhuis 0(n) class of models, with the 
weights for the vertex- interacting self-avoiding walk model (9 = w/2). K is the step 
fugacity and r = exp(— (3e) is the Boltzmann weight associated with non-intersecting 
site interactions. 
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walk. The choice p = then corresponds to a walk in which straight segments are 
eliminated (the walk is obliged to turn through a right angle after each step). The 
partition function for the model is then given by 

Z= £ (Kpf'p-^, (9) 
walks 

where r = exp(— j3e), N s is the length of the walk (number of steps), Nj is the number 
of doubly occupied sites and iV c is the number of corners. 

2.3. The corner transfer matrix renormalisation group method 

Following Baxter [SB E21 EB], the partition function of a two-dimensional lattice model 
may be written in terms of the product of four matrices representing the four quarters of 
the lattice, see Figure El The inputs and outputs of the matrices are the configurations 
at the seams of the four quarters. These matrices are known as corner transfer matrices. 
In general the four matrices are different, but may often be related by lattice symmetries. 
For our model the four matrices are the same up to a complex conjugation operation. 
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Figure 2. A schematic representation of the two-dimensional lattice split into four 
quarters, represented by the four corner matrices A, B, C and D. 




In general it is difficult to explicitly calculate these matrices for systems with a 
large number of sites. This is where the CTMRG method comes in^HJGI]; the matrices 
for larger lattices are calculated iteratively from smaller lattices. This is done as follows. 
A prototype system consisting of a small number of spins, four in the case of a vertex 
model, is set up exactly, see Figure El At each iteration the system is enlarged by adding 
spins, and this enlarged system is then projected back onto the prototype system in some 
optimal way, so as to minimise the loss of information. The spins of the prototype system 
have m states, where m determines the amount of information which may be carried 
forwards at each iteration. 

In addition to the vertex function W(a{, (jj, a^, ai), we also introduce the vertex 
functions P/v(ci, £j, £jfc) an d Cjv(£i, where the are the m-state spins of the prototype 
system. The matrix Cn corresponds to the estimate of the corner transfer matrix after 
N iterations. The partition function for a 2N x 2N lattice is given by 

Z N = Tr(C N C* N )\ (10) 

where the * denotes complex conjugation. The lattice is enlarged by introducing 
additional vertices, as shown in Figure 0] A corner matrix for the enlarged lattice 
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Figure 3. Figure (A) shows the initial vertex lattice. The 3-state spins defined in the 
vertex model are represented by circles. The black spins are summed over, taking into 
account the boundary conditions (fixed or free), to give the prototype system used in 
the CTMRG method, shown in (B). The squares represent the m-state spins used in 
the CTMRG method. Cn is an m x m corner transfer matrix representing one quarter 
of the lattice (see text). The asterisk represents complex conjugation. 




is given by 

C' N+1 = W*P N *P* N *C N , (11) 

where • indicates that the appropriate contractions over indices are performed. In a 
similar way we introduce an enlarged vertex function 

P' N+1 = W*P N . (12) 

By considering carefully the symmetry of the vertex function W it may be shown 
that the partition function for the enlarged lattice 2(N + 1) x 2(iV + 1) is given by 

Z'n+i = Tr(C^ +1 C^ +1 ) 2 . (13) 

We now wish to project the enlarged model, defined through C' N+1 and -Pjv+n back 
onto the smaller prototype system, giving Pn+i and CV+i- 
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Figure 4. The matrix C' N+1 is constructed by adding the vertex functions Pn and 
W to Cn as shown. The asterisk represents complex conjugation. Circles represent 
the 3-state spin variables defined by the original model, and the squares represent the 
m-state spins defining the CTMRG prototype system. The black spins are summed 
over, leaving the white spins as the indices of C' N+1 . C' N+1 is a 3m x 3m matrix which 
must then be projected in an optimal way to give the new (m x m) corner transfer 
matrix Cjv+i- 




Since the configuration space is smaller after projection, some information is lost. 
The trick is to minimise this information loss. This is done by performing a change of 
basis in some optimal manner. The partition function Z' N+1 is invariant by any change 
of basis, but the partial trace, giving Z^+i after the change of basis, is not. For a model 
with real weights it is clear that the optimal choice will minimise e = \Z' N+1 — Z N+1 \. 
This is done by diagonalising the matrix z = (C' N+1 C'^ +1 ) 2 and keeping the basis vectors 
corresponding to the largest eigenvalues. In a model with real weights z is just the 
density matrix, up to a constant. In the case considered here the eigenvalues of z 
are complex, coming in conjugate pairs. Minimising e corresponds to taking the basis 
vectors for the eigenvalues with the largest real parts, since the imaginary parts cancel. 
This may well involve dropping terms contributing to Z' N+l with larger moduli than 
other terms which are kept. It is not clear how this effect propagates with the number 
of iterations. Another choice is to keep the terms contributing to Z' N+1 with the largest 
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absolute values. This is found to give better results in practice. To complete the iteration 
scheme, we must also calculate Pjv+i- This is done by projecting P' N+1 onto Pn+i using 
the same basis vectors. As this iterative process is repeated, the partition function and 
other thermodynamic quantities are calculated for larger and larger lattices. The value 
of m defines the size of the configurational space kept from one iteration to the next. 
The larger the value of m, the better the approximation. 

The iteration scheme not only gives an approximation of the partition function for 
a square system, but also gives an approximation for the transfer matrix from which an 
infinite strip may be constructed. The transfer matrix is given directly in terms of the 
vertex functions Pjv and W as 

T = P N »W »W •P* N (14) 

for even lattice widths (width L = 2N + 2) and 

T = P N *W*P* N (15) 

for odd lattice widths (width L = 2N + 1). Free energies and correlation lengths for 
strips may then also be calculated, and the arsenal of finite-size scaling methods available 
for strips may be used. CTMRG is most efficient, however, when it is used to calculate 
directly one-point functions, such as the density, and this is how we use it in this article. 

The density is most simply calculated by selecting one of the central lattice bonds. 
The density p is then given by p = (H), where a is the value of the spin on the selected 
bond. In terms of the corner transfer matrices, 

_ Tr ( A C' N C%C' N C%) . . 

±T{Ls N Ls N U N U N ) 

where A is a matrix whose elements take the value \ a\ and L = 2N is the linear dimension 
of the lattice. 

In spin models the one-point functions are often calculated by the infinite lattice 
met hod 124*]. where the CTMRG method is iterated for a fixed value of m until the 
quantity of interest (magnetisation or energy per spin) no longer changes. This occurs 
because the number of eigenstates kept at each iteration is no longer sufficient to 
propagate the two point correlations throughout the system. The limiting value of 
the one point correlation function may then be studied as a function of m. In such a 
manner the magnetisation or energy per site may be estimated for the infinite system 
directly. This method relies on the rather nice way in which CTMRG breaks down for 
spin models. Unfortunately the complex nature of the vertex models at the heart of our 
calculations gives rise to more complicated behaviour as m is increased. Figure El shows 
the density of the self- avoiding walk model (p = 1, r = and n = 0, or equivalently 
9 = 7r/2) as a function of L for different values of m. It is clearly seen that there is no 
limiting value for the density as L is increased. Once the value of m is insufficient, p has 
a chaotic oscillatory behaviour with L. Instead of fixing m and varying L, we choose 
to do the opposite. Taking the density as an example, we define as the density 

for a system of linear dimension L for a given value of m. The limiting value of 
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Figure 5. Variation of the density p for the self-avoiding walk model (p = 1, r = 0) 
calculated using CTMRG as a function of the linear size L = 2N of the lattice for 
different values of m. 
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is calculated by increasing the value of m for fixed L until the required accuracy for 
Pl is reached. This is then studied as a function of L using finite-size scaling methods, 
described below. 

3. Results 

3.1. The vertex-interacting model with p = 

When p = the model is expected to have a critical behaviour equivalent to the standard 
O-point modeljlH|; for small enough r the average length of the walk diverges as a 
critical value of the step fugacity, K c , is reached. At K c the density is still zero and 
the gyration radius, Rq, (or any other characteristic measure of the size of the walk) 
diverges as K — > K~ following a power law defining the correlation length exponent v: 

£~(K C -K)- U . (17) 

For r smaller than some tricritical value, r tc , the walk is in the same universality class as 
the usual self-avoiding walk, and v = 3/4. This behaviour terminates at the tricritical 
point (K tc ,T tc ), where v = 4/7 as for the standard G-point. For values of r > r tc the 
length of the walk diverges discontinuously at a first order transition. The phase diagram 
of the p = vertex- interacting model is qualitatively the same as for the point [16J, 
and is shown schematically in Figure El Unlike for the standard 0-point model, the 
location of the tricritical point is known exactly to be (K tc = 1/2, r ic = 2)[16j. 

The idea is now to study this model using the CTMRG method outlined. As 
mentioned in sect ion l2~3*| CTMRG gives as a byproduct an approximation of the transfer 
matrix for a strip of width L through (fH^ or (|T3|). depending on the parity of L, 
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Figure 6. The schematic phase diagram for the vertex-interacting model with p = 0. 
For p = the vertex-interacting model has the same schematic phase diagram as 
the 0-point model, and the same critical behaviour. The solid line corresponds to a 
critical phase transition in the self-avoiding walk universality class. The dashed line 
corresponds to a first-order transition. The two types of phase transition are separated 
by a tricritical point (K tc — 1/2, r tc — 2). 
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from which the phase diagram and critical exponent u could be calculated using the 
phenomenological renormalisation group due to Nightingale [20 . This turns out, at least 
for the case at hand, not to be the most efficient way of approaching the problem, due 
to the additional computational effort required to calculate the relevant eigenvalues of 
T. The CTMRG method is far more efficient for the calculation of one-point functions, 
such as the bond density, which may be calculated in the middle of the system. With 
this in mind we use a phenomenological renormalisation scheme based directly on the 
scaling behaviour of density. 

The scaling ansatz states that the singular part of the free energy is a homogeneous 
function of its arguments (HI). This leads to the following finite size scaling expression: 

f s (K,L) = L- d f(\K-K c \L 1 /»)i (18) 

where d is the spatial dimension, here d = 2. Taking the first derivative, 

p s (K, L) = L^- 2 p(\K - K C \L^), (19) 
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In general 

p(K,L)= Poo (K)+p s {K,L), (20) 

and so it is necessary to know the value of p^K) in order to exploit fully the 
scaling behaviour. The density of a finite walk on an infinite lattice is zero, and so 
Poo{K < K c ) = 0. This enables the setting up of a phenomenological renormalisation 
group method using the function 

log (p s (K,L)/p s (K,L')) 
<-Pl,l'(K) = • (21) 

Using ()19j) we find that <£l,l' = 1/u — 2 independently of L and V . Naturally there are 
additional finite size corrections which should be taken into account, but the conclusion 
is that if the function (pL^>(K) is plotted for various values of L and V then it will 
converge to a fixed point given by (p(K c ) = 1/u — 2. In what follows we have set 
L' = L — 2 and looked for solutions of the equation 

Vl,l-2(K^) = ^, L ^{K L C ). (22) 

If such solutions exist then K c = lfm^oo and v = lfm^oo 1/(2 + ^l,l-2{K^)). 

Preliminary results using CTMRG were recently reported in reference [EJ; notably 
the density profiles and values for K c and v are given for various values of r with p = 0. 
When p = 0, the vertex-interacting B point is known to occur exactly at K tc = 1/2 and 
r tc = 2. In Figure U\ we present results using this scheme where we have fixed p = 
and t = r tc = 2. The results converge rather quickly to the exact values, K c = 1/2 
and z/0 = 4/7. Whilst here the tri critical values of K, t and v are known exactly, this 
is in general not the case. It is interesting to know how well the method works with 
no external input. As we pass through r tc the value of v changes from t'sAw — 3/4 for 
t < r tc to vq = 4/7 for r = r tc and on to v — 1/2 for r > r tc . This last value is simply 
a reflection of the fact that the density is nonzero, and so Rg ~ N 1 ^ 2 . Clearly the finite 
size estimates, vl{j\ will be continuous functions of r, tending to the step function 
only in the infinite lattice limit. If the different vl{ t ) cross, then these intersections 
necessarily tend, in the limit of an infinite system, to the correct values: v tc = ve and 
r tc . This is a standard method used to determine numerically the location, for example, 
of the point in the standard model (see for example [3 El EI])- The estimates derived 
from this method are shown in Figure El plotted as a function of 1/L. The estimates 
converge convincingly to the exact values r tc = 2, K tc = 1/2 and v tc = 4/7. 



3.2. The vertex-interacting model with p = 1 

What is interesting about the vertex-interacting model is that the phase diagram is 
very different when p ^ 0[35j. Naively one would expect the p = 1 model to also 
be equivalent to the 0-point model; the model consists of a walk which is essentially 
self-avoiding with short-ranged attractive interactions. This is in fact not the case; an 
additional transition arises in the dense walk phase and the nature of the transition 
at the (now) multicritical point changes. There are no exact results for the location 
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Figure 7. Estimates of the tricritical values K tc and v tc calculated for p = and 
r = Ttc — 2 using l|22|l . plotted as a function of the linear size L. There is a clear 
convergence to the expected exact values K tc = 1/2 and v tc = v$ — A/7 . 
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Figure 8. Estimates of the tricritical values (a) r tc , (b) v tc and (c) K tc calculated 
from the intersections of vl{j) for p — 0. There is a clear convergence to the exact 
values are r tc — 2, v tc — A/7 and K tc = 1/2 
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Figure 9. The density of the walk on the lattice for p = 1 and r = 1 as a function of 
K, clearly showing the existence of two phase transitions. 
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of the multicritical point (K mc , r mc ) as was the case for p = 0. There exist, however, 
exact predictions for the universality classes of the different phase transitions. The low- 
density phase transition is expected to be in the self-avoiding walk universality class 
(with v = 3/4) [30, 35j. The high-density critical line is expected to be in the Ising 
universality class {y = 1)[3U1|33]. The value of v mc at the multicritical point has been 
found from an exact result for the 19 vertex model with complex weights used in our 
calculations. It is predicted that v mc = 12/23^1- I n the remainder of this section we 
calculate the phase diagram and compare critical exponent estimates calculated from 
the CTMRG method with the known results. 

In figure 01 the density is plotted as a function of K for p = 1 and r — 1. The 
presence of two phase transitions may clearly be seen, the first between the zero-density 
phase and a dense walk phase, and the second between two dense walk phases. 

For the first transition, as K is increased, it is possible to use the scaling arguments 
given in section I3~T1 since the density in the thermodynamic limit is zero. For the second 
transition, however, the density is non-zero, and a priori unknown. To apply the scaling 
arguments it is necessary to estimate the density in the infinite lattice limit, poo(K). 
Alternatively it is possible to elimitate by using two lattice sizes. Along with (fT§j) 
and pU|) this gives a new finite-size scaling function 



which, as before, may be used to find estimates of the critical lines by looking for 




(23) 
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solutions of the equation 



<p L {K L c )=q> L _ 2 {K L c ). 



(24) 



In the limit L — > oo the solutions of this equation tend to the correct fixed points, such 
that (p(K%) -> l/v - 3. This again enables the calculation of the critical exponent v. 

The high-density transition line is expected to be in the Ising universality class, with 
an exponent v — 1. It is important to note that the exponent is used to represent both 
the critical exponent related to the correlation length and the inverse of the Haussdorf 
dimension, <i_f/fHB]- The Haussdorf dimension is the fractal dimension defined through 
the scaling of the density with the distance from the centre of mass of the walk. This 
leads to the scaling law: 



where Rq is the radius of gyration, or any other characteristic measure of the linear size 
of the walk. Since in the low density regime the gyration radius is proportional to the 
correlation length, then as the critical line is approached, and L —>■ oo, we may identify 
v = 1/dn- It is common practice to define v from (|25jl . whence z/=l/<i=l/2on 
the first-order transition line, where the fractal and space dimension become the same. 
The v associated with the high-density transition will necessarily not be defined in 
this geometrical fashion, but corresponds to the true thermodynamic critical exponent 
related to the correlation length. The estimate of K c and the associated estimates for 
v are given in Figure ITU1 for p = 1 and r = 1. The evolution of v is not monotonic, 
emphasising the need to go to large enough lattice sizes before having reliable estimates 
of the exponent. It is clear that v converges to the expected value v — 1, and we 
estimate the critical fugacity in this case to be K c = 1.326 ± 0.002. 

Figure UTT a) shows the phase diagram where both transition lines are calculated 
using the scaling function given in (J2~3~j) . In practice there are usually two solutions to 
(I23|) corresponding to the same transition; both converge towards the same value of 
K c as the size of the lattice is increased. However, one of these solutions evolves more 
slowly than the other as system size is increased, and this is taken as being the relevant 
solution for the calculation of v and K c . As r is increased these two solutions come 
closer (for fixed lattice size), and join at a certain value of r. Beyond this value of r 
there is nolonger a solution corresponding to the transition in question (there may be 
solutions corresponding to other transitions). This vanishing point moves to larger r as 
the lattice size is increased. The phase diagram shown in Figure HIT a) is calculated for 
system sizes up to L = 30. 

Figure HlT b) shows the estimates for the zero-density phase transition lines 
calculated using (|2T|) . As r approaches its multicritical value, the finite sized system 
feels the presence of the high-density phase transition, leading to a peak in the estimated 
transition line. This peak reduces in size rapidly as the system size is increased. 

Estimates of the location of the multi-critical point, {K mc ,T mc ), may again be 
found by looking for crossings in the finite-size estimates of u(t) calculated along the 
self-avoiding walk transition line using ()21)1 and ()22j) . These estimates are given in 



R G ~ N 1/dH , 



(25) 
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Figure 10. Estimates of K c and v corresponding to the high-density phase transition 
plotted as a function of 1/L for p = 1 and r = 1. The estimated value of 
K c = 1.326 ± 0.002 and convergence to the expected value of v = 1 is observed. 
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Figure 11. The phase diagram calculated using (a) l|23|) with L = 30 and (b) l|21|l . 
having fixed = 0, plotted for different values of L between 10 and 26. In (a) 
solutions corresponding to the low and high density phase transitions are shown, 
while in (b) only the low density transition was obtained. The circle in (a) shows 
the numerically estimated location of the multi-critical point (see text and Figure IT2"|) . 
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Figure 12. Estimates of (a) r mc , (b) v mc and (c) K mc calculated from the intersections 
of ^l(t) for the vertex- interacting model with p = 1 and fixing = 0. The estimated 
position of the multicritical point is then found to be K mc = 0.3408 ± 0.0002 and 
T~mc = 4.69 ± 0.01. The estimated value of v mc is incompatible with the exact value 
v mc = 12/23 w 0.5217 for reasons explained in the text. 
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figure [T2"l The location of the multicritical point is estimated as K mc = 0.3408 ± 0.0002, 
T mc = 4.69 ± 0.01. These estimates agree very well with the estimates given by Guo, 
Nienhuis and BTotejSHj, & ve z mc = iQn C T rac ~ 0.54, however the finite-size estimates 
of v^ir) cross at a value too low to be consistent with the exact value at the multicritical 
point (vl^oo ~ 0.51 compared with i/ exa ct = 12/23 = 0.5217- •■). This leads to the 
conclusion that either the scaling method is flawed or p^ ^ at the multicritical point. 
In Figure ED the density is plotted as a function of r along the critical (self-avoiding 
walk) line. At the multicritical point the density has a clear finite limit as the size is 
increased, indicating that the transition is first order as approached in this direction. As 
further evidence of the first order character of the transition when approached following 
the self-avoiding walk line we plot in Figure El the finite-size estimates of K c and p 
for different values of r near r mc . From Figure ^2 it is seen that for r close to r mc the 
system is influenced by the high-density transition. As the system size is increased, the 
estimates transit to the zone of influence of the zero-density transition, where they tend 
to their limiting values. This transit is clearly seen in Figure El for r = 4.3 < r mc . 
When t = 4.6 it is seen that the transit begins, but the sizes we could consider in 
a practical time scale were not large enough to see the end of the transition. It is 
however reasonable to suppose that for r = 4.6 the limiting density will be zero and 
occur as L — > oo. No such transit is seen for r = 4.7 > r mc , which displays a non-zero 
limiting density. The precise value of this limiting density is of course not reliable, since 
the location of the transition point, found assuming p^ = 0, will not be precise. The 
determination of the multicritical point is however expected to be correct since p^ = 
along the entire self-avoiding walk line, as must be the case since the fractal dimension 
d H = \jv < 2. 

There exist five branches of solutions for the Blote-Nienuis 0(n) model on the square 
lattice [3U|, through its mapping onto the 19 vertex model, for which the (multi) critical 
points are exactly known, as well as many of the critical exponents. The multicritical 
point here is thought to correspond to the same universality class as the n = solution 
belonging to the third branch. The location of this point is given exactly as[30j: 

z = K 2 t = {2- [1 -2sin(0/2)] [1 + 2 sin(^/2)] 2 } _1 ' 

K = -42sin(fl/2)cos(7r/4-#/4) > ^ 

pK = z[l + 2sin(0/2)] | 
6 = -it/4 

This gives the location of the exactly known multicritical point as K^ c = 0.446933 • • •, 
Pmc = 0.275899 • • • and r^ c = 2.630986 • • •. We checked the determination of this point 
using CTMRG fixing p = p* mc and r = r^ c and looking for the limiting value of Kl. This 
is shown in Figure El and is clearly consistent with the exact known value. The density 
was then calculated firstly at the finite-size estimate of the critical point and then fixing 
simultaneously the values of K, p and r to the values given in (J2H|) (see Figure ITB^) . In 
both cases a conservative estimate gives = 0.43 ± 0.01. These results indicate that 
the "multicritical" point is fairly exotic in that it corresponds generically to a first-order 
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Figure 13. The density plotted as a function of r calculated along the critical self- 
avoiding walk line K c {t). The density appears to be developing a jump at around r mc 
as the system size, L, is increased. 
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transition but with a critical behaviour if approached in a specific direction. We have 
tried to recover the critical exponent 12/23 by using the scaling behaviour of the density 
at K^ c , t^ c and p* mc . Whilst the results are not inconsistent with v = 12/23, they are 
not of very good quality and are also not inconsistent with a value v « 0.54 as originally 
found by Blote and Nienhuis using transfer matrices jHUj. Whilst in principle there is no 
limit to either the precision or lattice sizes obtainable using CTMRG, there is a practical 
limitation in terms of computing time and resources. 

4. Conclusion 

In this paper we have presented a version of the corner transfer matrix renormalisation 
group method suitable for use with self-avoiding walk type models and related 0(n) 
symmetric models. We have applied it to the vertex-interacting self-avoiding walk 
model. The method was shown to reproduce the known results, and has given new 
insight into the critical nature of the model when p ^ 0. 

Self-avoiding walk models are used as models for polymers in dilute solution, and 
it is known that the lattice #-point model provides a good description of what happens 
for real homopolymers in solution^]. The standard textbook reason for not taking into 
account the rigidity of the polymer is that the rigidity simply changes the persistence 
length, or in terms of the lattice model the number of real monomers per step of the walk. 
The critical behaviour is unaffected in the infinite walk limit. This should be true as 
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Figure 14. (a) Shows the estimates of K c and (b) estimates of the density plotted 
as a function of 1/L for three values of t: t = 4.3 and r = 4.6 which are below 
T mc and t = 4.7 which is above r mc . When r < r mc the graphs show a transit from 
being influenced by the high-density transition line (small L) to being influenced by 
the zero-density transition line (large L). Whilst the transit is not complete for the 
sizes considered when r = 4.6, it is clear that = in this case whilst ^ when 
t = 4.7. The solid line is given as a guide to the eye. 
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Figure 15. Estimates of K c calculated using Q21|l with p = p* mc and r = r^ c fixed at 
the values given by the exact solution (|26(l . The estimates converge well to the exact 
value K* nc = 0.446933 • • • (shown by the solid line) as L — > oo. 
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Figure 16. The density plotted for (a) p = p* mc and r = r^ c given by the exact 
solution H26J1 whilst K mc is calculated from (|21|l (filled circles) and (b) p, r and K 
fixed by (|26|l (squares). The density does not become zero in the limit L — > oo. A 
conservative estimate of the limiting density would be 0.43 ± 0.01. 
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long as the fractal dimension dn is smaller than the spatial dimension (i.e. d# < 2 here). 
However, when djj = d it may be expected that the polymer feels the presence of the 
lattice. This will then give rise to competition between the rigidity and the underlying 
lattice structure, reminiscent of frustration in spin models (indeed the critical behaviour 
of the Nienhuis 0(n) model on the square lattice has been related to the fully-frustrated 
XY model[2n])- If such a claim is true, then the critical behaviour observed will be 
universal (in the sense of not depending on the lattice) and independent of rigidity if 
dn < d, and may be non-universal (in the same restricted sense) and may depend on 
the rigidity if djj = d. One could therefore argue that the simple observation that 
the critical behaviour of the collapse transition changes as p is varied implies that the 
density is non-zero at the collapse transition when p ^ 0. In interpreting this transition 
as due to frustration, it should be noted that it does not exist for the 0(n) model on 
the hexagonal lattice [2fi| 127]. 

There are other models which have qualitatively similar phase diagrams. Most 
well known is the self-avoiding trail model [3*7] . A trail is defined to be self-avoiding on 
the bonds, but is allowed to visit sites more than once, and in particular is allowed to 
intersect itself. The attractive interaction is again associated with doubly visited sites 
(intersection or collision). It has been conjectured that the collapse transition for trails 
occurs exactly at r = 3 with an exponent v = 1/2; it is also conjectured that limiting 
density is zero. This conjecture resulted from a Monte-Carlo enumeration of extremely 
long walks of the kinetic trail model, believed to be equivalent to the trail model at the 
collapse transition [HH] • As discussed above, in this case what is being calculated is in 
fact the fractal dimension djj- The relation v = 1/dn is only valid when the critical 
length may be measured using the gyration radius. As argued above, the validity of 
this identification may be questioned when du = d. It would be of interest, then, to 
check by other means the possibility that the criticality is in fact not of the same type 
as observed here for p ^ 0. This question is currently under investigation. 

Another model with a qualitatively similar phase diagram is the Hydrogen-bonding 
model (T7J. This model includes interactions between nearest-neighbour sites, which may 
be taken into account in the CTMRG method by mapping the walk model onto a 4-state 
vertex model with complex weights. There are, however, important differences, notably 
the high-r dense- walk phase has p = 1 throughout, whereas here p ^ 1 except in the 
limit t — > oo. To what extent this may affect the nature of the phase transitions is not 
clear; early results indicate that perhaps the dense-walk phase transition is not in the 
Ising universality class [39] . 
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